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Abstract 

We investigate a lattice Nambu-Jona-Lasinio model both by the 
Monte Carlo method and Schwinger-Dyson equations. A compari- 
son allows the discussion of finite size effects and the extrapolation 
to infinite volume. We pay special attention to the identification of 
particles and resonances. This enables us to discuss renormalisation 
group flows in the neighbourhood of the critical coupling where the 
chiral symmetry breaking phase transition takes place. In no region of 
the bare parameter space do we find renormalisability for the model. 



1 Introduction 

The study of four fermion interactions has a long history. Fermi originally 
introduced the idea as a theory for weak interactions |1|. Later Nambu 
and Jona-Lasinio proposed their model (in the following called the NJL 
model) to describe spontaneous chiral symmetry breaking and dynamical 
mass generation. This was one of the first works to note the association 
of a massless Nambu- Goldstone boson with the spontaneous breaking of a 
continuous symmetry 0. 

Over the years there have been efforts to investigate the vacuum struc- 
ture of strongly interacting four fermion theories 0]. Recently there has been 
renewed interest in this subject in connection with technicolour and the top- 
mode standard model ^, |^. For these applications one is mainly interested 
in the physics of the broken phase where the fermion mass is dynamically 
generated. Additional interest was generated by the discovery that a general- 
ized NJL model is equivalent to the Higgs- Yukawa model 0, ^ and possibly 
to SU{N) gauge theory |^. Furthermore there have been various investiga- 



tions of the gauged NJL model |]T0| which show a line of second order phase 
transitions linking the pure NJL model and QED. This raises the question 
of whether the two models fall into the same universality class. Our investi- 
gation of QED [|ll[] has shown that in the strong coupling regime (including 
the critical region) it is impossible to vary the cutoff without changing the 
physics, while in the weakly coupled region this is possible. We have called 
such behaviour weak renormalisability. 

In this paper we are mainly interested in the model's phase transition and 
critical behaviour. Since this falls into the large coupling region, nonpertur- 
bative methods are required. We use the lattice Monte Carlo simulation 
and the semi-analytic Schwinger-Dyson method [0, |13|. The Monte Carlo 
method has the advantage of being exact but can only be applied on finite 
lattice volumes. The Schwinger-Dyson approach on the other hand can be 
extended to infinite volumes but suffers from having an uncontrolled approx- 
imation. The synergy of both methods allows us to go much further than 
with either method alone. 

We have looked at the chiral condensate in order to map out the phase 
diagram. We find a second order chiral phase transition implying infinite 
correlation length which in principle allows us to obtain a continuum theory. 
In the vicinity of the critical point we calculate the renormalised fermion 
mass (/ir) and the energy levels of the pseudoscalar (vr, Goldstone boson), 
scalar (a) and vector (p) fermion-antifermion composite states as well as the 
pion decay constant (/t^). 

A major problem is the identification of resonances and bound states on a 
finite lattice. For the vr and p we were able to identify the first excited energy 
levels by our Monte Carlo calculation. These compare favourably with the 
results from the Schwinger-Dyson equations. Thus we feel encouraged to take 
the thermodynamic limit of the Schwinger-Dyson equations. This allows the 
calculation of the spectral functions and thus the identification of resonances. 
With this identification at hand we construct flow diagrams of dimensionless 



ratios of physical quantities. We use them to investigate the renormahsation 
group flow in the vicinity of the critical point. 

The paper is organized as follows. Our lattice model is defined in section 
2.1. Details of the Monte Carlo method are given in sections 2.2, 3.1, 3.2 



and 4.1 (see also fl^)- The Schwinger-Dyson equations are discussed in 
sections 2.3, 3.2, 3.3, 4.2, 4.3 and the appendices. Our main results of both 
methods are given in section 5 on meson spectroscopy and in section 6 for the 
renormahsation group flows. Finally in section 7 we present our conclusions. 



2 The Model 

Because we are interested in chiral symmetry properties we use staggered 
fermions throughout. 

2.1 The Action 

The NJL model is defined by the continuum action 



S = d X \ i){x)'^^d^^{x) - qq 



(i){x)i){x)) - (V^(x)75^(x)) 



(2.1) 



which is chirally invariant. For staggered fermions we take the action 

^ X,fl X 

- fl'oZlx(a;)x(a;)x(a; + A)x(a; + A), (2-2) 

where 

i]^[x) = (_i)-i+-+-.-i, r],{x) = 1. (2.3) 

The interaction is chosen as the simplest one that is invariant under the 
continuous chiral transformation 

X{x) -. e^"^(^')x(x), x{x) -^ e^"^(^)x(a:), (2.4) 

with 

e(x) = (_i)-i+-2+-3+x4_ (2.5) 

The mass term acts as a chiral symmetry breaking external source and as an 
infrared regulator. 

2.2 Monte Carlo Simulation 

In order to use the Monte Carlo method we first have to integrate out the 
fermion fields. Therefore we have to rewrite the interaction in a bilinear 
form. This is done by introducing an auxiliary field Ofj,{x) G [0,27r) |l^, |T6|: 

-^int = V9~oT.V,ix) [x{x)e''^^'^^xix + fi)-xix + i:i)e-''^^'hix)] . (2.6) 

X,fl 

The Hybrid Monte Carlo method was employed to generate the configura- 
tions. Lattices of size V = L^ x Lt = Q^ x 12, 8^, 8^ x 16, 8^ x 32 and 12^ 
were used. Boundary conditions were periodic in the spatial directions and 
anti-periodic in the time direction. We investigated the model for niQ = 0.01 
to 0.09 and Qq = 0.21 to 0.32. For each parameter set 1000 trajectories were 
generated for measurements. The trajectory length was always ~ 0.8. The 
number of time steps was adjusted such that acceptance rates between 70 
and 80% were achieved. Propagators were calculated with the conjugate gra- 
dient method. As convergence criterion we required that the norm of the rest 
vector was smaller than 10~^. The chiral condensate was calculated with a 
stochastic estimator (see, e. g., ||T7|). 



2.3 Schwinger-Dyson Equations 



The model ( |2.2|) was also studied by solving truncated Schwinger-Dyson 
equations up to order go and Qq on the lattice. To make a direct comparison 
with Monte Carlo results possible we use the same action with the same cut- 
off (for another application of Schwinger-Dyson equations to a lattice theory 
see H). 

In order to write the equations in this section in a compact form it is 
useful to give the lattice action (p.2|) with a slightly different notation: 



^ = E X.d^yX'^ + mox.x" - goXxX^'cl^^'XyX'^' ■ (2.7) 

The fields at site x are Xx and x^^ ^^^ the Einstein summation convention 
is to be used whenever a site index occurs twice. The meaning of the kinetic 
matrices d^ and of the coupling tensor c^,y should be clear by comparison 
with (^I^) : d^^ y = +\ri^{x) ii y = x + fi, d^ y = —\rj^{x) ii y = x — ft and 
zero otherwise; c^,y, = 1 ii x = x', y = y' and x and y are neighbours, and 
is zero otherwise. It is convenient to give the coupling tensor four indices so 
that we can use the Einstein summation convention and to remind us that 
more general four-fermi coupling terms are possible. 

The lattice Feynman rules are that the bare staggered fermion propagator 
is 

G''^^^=(moS;; + Y.d^A , (2.8) 

and the vertex is 

- 9oc:,'y, . (2.9) 

When the interaction is turned on the bare fermion propagator is replaced 
by the full fermion propagator 



Gf"ii= M^^ + E^^^ • (2.10) 



The new matrices My and D'^ y are generalisations of the matrices mo6y and 
d'^ y appearing in (|2.8|) and have the same symmetries, e. g.. My is symmetric 
under reflections and links sites with even separations in violation of chiral 
symmetry. On the other hand D^ y is odd in the /i direction and even in 
the other directions and only links sites with odd separations and so respects 
chiral symmetry. M is invariant under translations, D'^ acquires the same 
factors of ±1 as c?'^. 

In fig. we show graphically the Schwinger-Dyson equations for the 
fermion propagator. It is not possible to solve the full Schwinger-Dyson 
equations. They have to be truncated at some stage. The crudest approxi- 
mation is to keep only term (a) which gives the gap equation, namely 

^ y ~ " ?; ' 
M^ = N5^, (2.11) 
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Figure 1: Pictorial representation of the Schwinger-Dyson equation for the 
ferniion propagator. 



where N satisfies 



N 



mo 



1 



^077 E 



N 



V 



(2.12) 



The /c-sum runs over all momenta consistent with the boundary conditions. 
The fermion mass /ir is given by the position of the pole and is /iR = 
sinh"^ A^. 

The next approximation is to keep terms (a) and (b). Term (a) leads 
simply to a renormalisation of the mass and term (b) leads to a renormalisa- 
tion of the kinetic energy. We shall call these the order go Schwinger-Dyson 
equations because the error is of order Qq. Explicitly they read 



y 

My 






(2.13) 



where iV and F satisfy 



N 



"^o + S^ottE 



N 



V^ N^ + Y.uFisin'k, 



2i, ' 

V 



i+2<,„ii:- 



F^sin /c^ 



E.F^in'A;, 



The pole position giving the fermion mass is 



m 



sinh — . 
Ft 



(2.14) 



(2.15) 



In infinite volumes (when the sum is replaced by an integral) F needs no 
direction index, as all the F^'s are identical. However in a finite system 
with different lattice lengths and boundary conditions in the space and time 
directions an index is needed. 

The order g^ equations also include term (c), but with the full vertex 
replaced by the bare vertex, an approximation which introduces an error of 
O^Qq). The resulting equations (given in configuration space) are 

+ doCyv^ u~90Cvy^ u 

2 xu/^iull r/^fuU Z/^fuU v sw 

yQ^vr s u w zy 

I ^2 xu/^iuU r /^fuU Z/^fuU v sw fry -i a\ 

M and D can be projected out by using their different behaviours under 
reflection. For the flrst time the unknown matrices M and D are not sim- 
ply renormalised versions of the bare matrices, but extend further in space. 
Mass and wave- function normalisation become momentum dependent. The 
matrices M^ and D^ are stored in position space. These individual compo- 
nents can be thought of as the Fourier coefficients for the momentum space 
fermion propagator. 

The equations were solved numerically by simple iterations. Typically 
about 7 iterations were needed in the order qq case to obtain convergence 
and 20-40 iterations in the order g^ case. 



3 Chiral Condensate and Fermion Mass 

3.1 Phase Diagram 

To find the phase diagram we have computed the chiral condensate (xx) 
also referred to as a. Monte Carlo results are collected in tables |l| and ^ 
To show that we are working in the vicinity of the critical point we look at 
a Fisher plot ||l9l of our data, i. e., we have plotted {xxY against mo/ {xx) 
in fig. ^ The lines connect points of constant g^. If the theory was mean- 
field like, these would be straight parallel lines [|19|. We see that this is 



approximately the case for our data. For large qq the extrapolated lines end 
on the y-axis giving a finite value for (xx) at ttiq = 0, while for small Qq the 
extrapolated lines end on the x-axis giving zero for (xx) at rriQ = 0. The 
line that corresponds to the critical coupling would end at the origin as is 
indicated by the dashed lines in fig. ^ On the 8^ x 16 lattice this gives an 
estimate of about 0.25 < Qc < 0.26 while on the 12^ lattice Qc is somewhat 
larger than 0.27. We see from this that there are finite size effects in the data 
as can also be inferred from a cursory look at tables |I] and |^. We shall in the 
following call the region go < g^ the symmetric phase and the region go > gc 
the broken phase. 

In fig. ^ we compare our Monte Carlo data with results from the gap equa- 
tion and our two Schwinger-Dyson approximations. (The chiral condensate 
is the value of the fermion propagator at zero distance: a = V~^ J2x G*^"" x-) 
There is a large discrepancy between the numerical data and the gap equa- 
tion. However the agreement improves considerably as higher orders are 
included leading to a maximum discrepancy of about 10% with our Monte 
Carlo data. The equations also reproduce the large finite size effects. 

Applying the Schwinger-Dyson equations for infinite volume and zero bare 
mass reveals that there is a chiral phase transition and not just a crossover. 
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Figure 2: Fisher plots of Monte Carlo data for the chiral condensate a on 
(a) 8^ X 16 and (b) 12^ lattices. 
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Figure 3: Comparison of Monte Carlo data (8^ x 16 lattice, rriQ = 0.02J for the 
chiral condensate with results from the gap equation (dotted line), the order 
go (dashed line) and the order g^ (solid line) Schwinger- Dyson equations. 



In the gap equation this occurs at 



1 



— ^0.201700, 
8y 



(3.1) 
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8^ 12^ oo 


0.21 


0.01 


0.0294(1) 0.0284(1) 0.0275(3) 


0.23 


0.01 


0.0419(3) 0.0399(1) 0.0387(5) 


0.25 


0.01 
0.02 


0.0632(8) 0.0616(4) 0.0608(11) 
0.1044(10) 0.1005(3) 0.0992(8) 


0.27 


0.01 
0.02 


0.1039(12) 0.0999(10) 0.0986(18) 
0.1421(13) 0.1358(5) 0.1343(10) 


0.28 


0.01 


0.1283(11) 0.1242(9) 0.1230(15) 


0.30 


0.01 
0.02 


0.1675(20) 0.1683(6) 0.1685(10) 
0.1929(12) 0.1896(4) 0.1892(6) 


0.32 


0.01 


0.2021(11) 0.2042(7) 0.2045(9) 



Table 1: Monte Carlo results for the chiral condensates on lattices with 
Lg : Lt = 1 : 1. The third column gives the extrapolation to infinite volume 
using eq. p.2-^ j. 



ixx) 


9o 


rriQ 


6^ X 12 


83 X 16 


oo 


0.21 


0.02 


0.1100(7) 


0.0784(2) 


0.0559(18) 
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0.1458(8) 


0.1170(4) 


0.0963(19) 




0.09 
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0.1794(2) 
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0.23 


0.02 


0.1354(10) 
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0.1182(17) 
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0.1915(7) 


0.25 
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0.0991(7) 


0.0720(43) 
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0.1604(8) 


0.1271(5) 


0.1034(25) 
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0.1890(6) 


0.1627(5) 


0.1471(17) 
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0.2339(6) 


0.2185(3) 


0.2124(8) 


0.26 
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0.1721(7) 


0.1412(6) 


0.1205(27) 
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0.2014(12) 


0.1757(4) 


0.1618(18) 
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0.2278(3) 


0.2224(7) 
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0.01 


0.1592(13) 


0.1324(8) 


0.1139(38) 
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0.1853(9) 


0.1555(5) 


0.1369(24) 
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0.2096(8) 


0.1878(5) 


0.1767(16) 
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0.2483(9) 


0.2364(3) 


0.2325(8) 


0.28 


0.02 


0.1962(9) 


0.1700(5) 


0.1542(24) 




0.04 


0.2197(9) 


0.2004(4) 


0.1915(14) 
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0.2560(6) 


0.2452(3) 


0.2419(7) 


0.29 


0.02 


0.2068(12) 


0.1846(7) 
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0.2305(8) 


0.2123(3) 
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0.2621(5) 
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0.1996(6) 
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0.2235(5) 


0.2180(13) 
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0.2698(6) 


0.2605(3) 


0.2580(6) 


0.32 


0.02 


0.2356(10) 


0.2225(6) 


0.2172(17) 
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0.2531(9) 


0.2433(5) 


0.2404(11) 




0.09 


0.2798(6) 


0.2728(4) 


0.2711(7) 



Table 2: Monte Carlo results for the chiral condensates on lattices with 
Ls : Lt = 1 : 2. The third column gives the extrapolation to infinite volume 



using eq. ( |3.2^ ). 



where 



y 



2- d^k 



(27r)4 EnSin'k, 



and in the order g^ case at 

9c = 



roc 

/ da e-^"I^{a/2) ^0.619734, (3.2) 

Jo 



32y 



ilQy-1] 



0.249483 . 



(3.3) 
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Figure 4: Monte Carlo data for the chiral condensate and curves resulting 
from the phenomenological fit eq. (\3.^ ). The dashed line indicates the extrap- 
olation to infinite volume and zero bare mass. Symbol shape denotes bare 
mass (diamonds 0.09, triangles 0.04, squares 0.02 and circles O.Olj. Colour 
denotes lattices size (white 8^ x 16 and black 12^j. 



In the order Qq case we can give no formula but the number is Qc ~ 0.263. 

The semi-quantitative success of these equations motivated us to make 
a phenomenological fit to the Monte Carlo data based on a modified gap 
equation 

C5^ M 



a 



v\ M^ 



J2u sin^ k^ 



where a is the chiral condensate and 



M = (ci + 025-0 + cafi'o)^- + C4mo . 



(3.4) 



(3.5) 



The values of the parameters were found to be: Ci = 0.149(1), C2 = 6.37(1), 
cg = -3.07(1), C4 = 0.938(1) and C5 = 0.959(1). The resuhs of this fit 
are plotted in fig. | which includes data from different lattice sizes showing 
that the fit equation also accommodates quite nicely the finite size effects. 
Extrapolating the phenomenological equation to infinite volume and zero 
mass gives Qc ~ 0.278. Comparing this value with the above results gives a 
shift of about 0.015 relative to the order g^ Schwinger-Dyson equations. 

The Schwinger-Dyson equations as well as the fit equation in the infinite 
volume limit lead to the following scaling laws sufficiently close to the critical 
point: 

a'^\n{ao/a) (xgo- gc (3.6) 
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Figure 5: Monte Carlo data for the chiral condensate extrapolated to in- 
finite volume and curves resulting from the mean-field fit with logarithmic 
corrections eq. I\3.1J\) . The dashed line indicates the extrapolation to zero 
bare mass. Symbol shape denotes bare mass (diamonds 0.09, triangles 0.04, 
squares 0.02 and circles O.Olj. 



on the critical line mg = and 

cr'^ln(cro/cr) oc mo 



(3.7) 



at Qo = Qc (co is a constant of order 1). This is mean-field behaviour with 
logarithmic corrections. Although at this order the power of the logarithm is 
+1 it is possible that the inclusion of higher order Feynman diagrams could 
change this power. 

For a more direct determination of the critical exponents and logarithmic 
corrections we extrapolate the data to the infinite volume first and then fit 
the result with a simple effective potential. The extrapolation of the chiral 
condensate to the infinite volume is done in section 3.3, and the results are 
stated in tables 1 and 2. Wherever we have two extrapolations, we take the 
one coming from the larger volume. We have performed two fits. In the first 
fit we assume an equation of state of the form 



mo = {5- - + l)Ka 



5-- 



+ {5 + l)Ca^ 



with 



K 

c 



Co + Ci(^o- 



(3.8) 



(3.9) 
(3.10) 



12 



thus leaving the critical exponents as free parameters. We find ki = —3.1(1), 
Co = 2.6(2), Ci = 8(1), /3 = 0.42(1), 5 = 3.4(1) and g^ = 0.280(1). The 
exponents f3 and 6 lie rather close to their mean field values (1/2 and 3 
respectively). In the second fit we fix f5 and 6 to their mean field values but 
allow for logarithmic corrections: 



tuq = 2k, 



Inkier 



:n + 4C 



In^ 



a 



(3.11) 



The result of the fit is ki = -3.6(3), Co = 2.1(3), Ci = 9(2), p = 0.36(11), 
q = 0.84(1) and gc = 0.282(2). The data and the fit are plotted in fig. g 
(The former fit leads to similar curves.) It is pleasing to note that all three 
fits give roughly the same value of gc- It is however possible that the powers 
p and q found from the fit eq. ( p.ll| ) have not yet reached their asymptotic 
values in the parameter region we have covered. 

To gain a first impression of the short distance behaviour of the model 
we have taken our data and in fig. ^ plotted the four fermion condensate, 
<74 = \J2x,fi{x{^)x{^)x{^ + f'')x{^ + P'))/V, against the chiral condensate. 174 
was also calculated via a stochastic estimator, but now involves two fermion 
matrix inversions. (The number of sets of random numbers used and hence 
inversions performed was also increased from one for the chiral condensate to 
~ 0(10%) of the number of molecular dynamics steps in the Hybrid Monte 
Carlo program.) The results are given in table ^. The points appear to lie 
roughly on one curve. We have compared the Monte Carlo data with the 
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Figure 6: Monte Carlo results for the four fermion condensate plotted against 
the chiral condensate. The curve results from eqs. ( |3.J2D and (\3.13i ). The 
symbols are the same as in fig. ^. 
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0.083(4) 


0.091(1) 


0.1058(6) 
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0.086(4) 


0.088(3) 


0.097(2) 


0.1073(5) 


0.28 


0.081(7) 




0.092(3) 


0.096(1) 


0.1105(4) 


0.29 






0.089(4) 


0.102(1) 


0.1133(6) 


0.30 


0.080(8) 


0.102(3) 


0.106(5) 


0.102(1) 


0.1156(5) 


0.32 


0.088(8) 




0.103(4) 


0.111(3) 


0.1198(6) 



Table 3: Monte Carlo results for the four fermion condensate a^. 



naive result obtained from free fermion propagators. This gives 



a4{mj 



E 



sm pi, 



-y- 

V p EAsin'^pA + m2_ 



am 



— [1 — mcr(m)] + a{mf' 
16 



(3.12) 



against 



cr(m) 



m 



E 



^ p EAsin Pa + 



m^ 



y^oo 



d'^p 



m 



^ (27r)4 ^;, sin^A + "^^ 



(3.13) 



where m is to be regarded as a parameter. The infinite volume result is also 
plotted in fig. |^. The data lie close to the upper branch (this corresponds to 
small /iR, while the lower branch corresponds to large /xr). This similarity 
between the data and the free fermion formula suggests that there is only 
little interaction between the fermions, the main effect of the interaction 
being the generation of a renormalised mass, the determination of which we 
shall now turn to. 



3.2 Fermion Propagator 

To determine the fermion mass the fermion propagator 

G{t)= Y. (X(x,t)x(0,0)) 



(3.14) 



x^ .xo.xo even 



is fitted by 
G{t) 



[I + e-'^R-^*) cosh/iR 



-Mr* _ C_]^yg-MR(it-t) 



(3.15) 
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Figure 7: Fermion propagators from Monte Carlo calculations on an 8^ x 32 
lattice at mo = 0.04, go = 0.23 (a) and go = 0.28 (b). Open (solid) symbols 
represent positive (negative) values. The lines connect values from fits using 
eq. ( \3.l5i) . The resulting fermion masses are (a) fj,R = 0.262(1) and (b) 
fXR = 0Al7{l). 



The parameter /ir is the renormahsed fermion mass and Z2 is the wave- 
function renormahsation constant. The results of the fits are given in tables ^ 
^ and ||. In fig. |^ we show two propagators on 8^ x 32 lattices which were 
simulated to check the masses obtained on lattices with smaller time extent. 
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Table 4: Monte Carlo results for the wave-function renormahsation constant 
Z2. 
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Figure 8: Monte Carlo results for the chiral condensate plotted versus the 
renormalised fermion mass. The dotted lines are results from the gap equa- 
tion. The solid lines are results from the order g^ Schwinger- Dyson equations. 
The symbols are the same as in fig. ^ 

We see that their behaviour is very well described by a single exponential. 

In fig. H we plot {xx) against /ir. In this plot we have also shown the 
predictions of the gap equation and of the order Qq Schwinger-Dyson equa- 
tions. We see that the data lie in a narrow band. The theoretical curves 
also display a small width because they still have weak dependences on bare 
mass and lattice size. The small spread tells us that the form of the fermion 
propagator is hardly changed by the interactions. 

3.3 Finite Size Effects 

Comparing results found on lattices with different volumes (see tables |T| 
and 0) makes it clear that we have large finite size effects in our data, and 
that it is important to study these in order to see the picture for an infi- 
nite lattice. An advantage we have over the usual situation in a Monte Carlo 
calculation is that we can use the Schwinger-Dyson equations to calculate re- 
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Table 5: Monte Carlo results for the renormalised fermion mass on lattices 



with L, 



volume using eq. (|3.25| ) 



1:1. The third column gives the extrapolation to infinite 



suits on large or even on infinite lattices, and so see how the thermodynamic 
limit is approached. 

To keep the calculations reasonably simple we will look at the lattice size 
dependence of the solutions of the order go Schwinger-Dyson equations ( p.l4| ), 
but we expect that the form of the formulae found will remain true in general. 

The equations ( |2.14| ) always involve the sum over allowed momenta of 
some function a{k) of the momentum. Let us call such a sum Ua- We need 
to compare this sum over momenta with the integral over all momenta which 
gives the infinite volume limit. We do this by using the Poisson resummation 
formula to rewrite our sum as an integral 



U, 



V 



k 

-' ^ k 

The sum of 6 functions can be expressed as a Fourier series, leading to 

Ua = J2i-^rMnLs,ntLt), (3.17) 



where A is the Fourier transform of a defined by 



A{x) 



{2ny- 



e''-''a{k) 



(3.18) 



(the (—1)"* arises because of the antiperiodic boundary conditions in time). 
The functions that we consider will always have Fourier transforms which 
decrease exponentially at large distance. Therefore in an infinite volume only 



17 



/iR 


9o 


mo 


63 X 12 


8^ X 16 


oo 


0.21 


0.02 


0.183(1^ 


0.133(0^ 


0.098(3) 




0.04 


0.250(1^ 


0.209(0^ 


0.179(2) 




0.09 


0.377(1^ 


0.345(0^ 


0.328(1) 


0.23 


0.02 


0.240(2^ 


0.181(0^ 


0.138(4) 




0.04 


0.311(1^ 


0.258(0^ 


0.222(3) 




0.09 


0.430(1^ 


0.399(0^ 


0.385(1) 


0.25 


0.01 


0.268(3^ 


0.185(1^ 


0.123(8) 




0.02 


0.305(2^ 


0.238(1^ 


0.191(4) 




0.04 


0.364(2^ 


0.317(1^ 


0.289(3) 




0.09 


0.483(1^ 


0.455(1^ 


0.445(2) 


0.26 


0.02 


0.334(2^ 


0.273(1^ 


0.232(6) 




0.04 


0.397(2^ 


0.352(1^ 


0.328(3) 




0.09 


0.512(1^ 


0.485(1^ 


0.475(2) 


0.27 


0.01 


0.329(3^ 


0.261(1^ 


0.213(8) 




0.02 


0.365(2^ 


0.306(1^ 


0.269(5) 




0.04 


0.430(2^ 


0.382(1^ 


0.358(3) 




0.09 


0.538(1^ 


0.515(1^ 


0.508(2) 


0.28 


0.02 


0.416(3^ 


0.345(1^ 


0.303(6) 




0.04 


0.462(2^ 


0.418(1^ 


0.398(3) 




0.09 


0.572(2^ 


0.548(1^ 


0.541(2) 


0.29 


0.02 


0.442(3^ 


0.385(1^ 


0.355(6) 




0.04 


0.492(2^ 


0.455(1^ 


0.440(2) 




0.09 


0.593(2^ 


0.574(1^ 


0.569(2) 


0.30 


0.01 


0.463(5^ 


0.387(2^ 


0.343(10) 




0.02 


0.476(3^ 


0.432(1^ 


0.412(5) 




0.04 


0.534(3^ 


0.492(2^ 


0.477(4) 




0.09 


0.630(2^ 


0.608(1^ 


0.602(2) 


0.32 


0.02 


0.563(5^ 


0.501(2^ 


0.475(7) 




0.04 


0.585(3^ 


0.562(1^ 


0.555(3) 




0.09 


0.681(2^ 


0.658(1^ 


0.653(2) 



Table 6: Monte Carlo results for the renornialised fermion mass on lattices 
with Ls : Lt = 1 : 2. The third column gives the extrapolation to infinite 



volume using eq. (|3.25| ). 



A{0) survives, while on a large lattice the finite size effects will be dominated 
by terms where one of the n^ is ±1. Physically speaking these represent 
the contributions of particles which have travelled "once around the lattice" . 
Terms where two of the n^'s are non-zero, or one of the n^'s is larger than 
1 die away as exponentials with higher exponent than the leading correction 
and very quickly become negligible. 

To see what the finite size corrections to the fermion propagator are we 
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Figure 9: Finite size scaling of the chiral condensate. Plotted are Schwin- 
ger-Dyson results on lattices with Lg > Q; Lg : Lt = 1 : 1 (solid symbols) 
and Lg : Lt = 1 : 2 (open symbols). The parameters are tjiq = 0.02 and 
^0 = 0.23,0.25,0.28. 



need to know the large distance behaviour of the Fourier transforms of 



N 



and 



N^ + EuF^s^^ K 



Ff, sin kf. 



(3.19) 



(3.20) 



i. e., the large distance behaviour of the fermion propagator. These limits 
can be found by the saddle point approximation. For definiteness we look at 
the case rii = 1. First integrate over ki which can be done exactly by looking 
at residues at the poles of the integrand: 
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Figure 10: Finite size scaling of the chiral condensate (Monte Carlo results). 
The lattice sizes are 6^, 8^ and 12^ (•) and 4^ x 8, Q^ x 12 and 8^ x 16 (o). 
The parameters are itlq = 0.02 and go = 0.25, 0.27, 0.30. 



d'^k 



Lseven / dk2dk'^dk/i 



exp{ikiLs 



N 



1 



N 



X exp 



— L, sinh 



(2^)' Jn^ + E.^1 F^ sin2 A;, JiV^ + Ff + E.^i F^ sin^ k. 



(3.21) 






For large Lg this integral is dominated by sin /c;^ ^ and is approximately 
Gaussian: 



^) F2F3F4 



oc (/Ur/L J 2 exp(-/iRL<,) . 



exp 



-L,sinh"'(iV/Fi 



(3.22) 
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(In the final expression /ir -C 1 has been taken.) The change in the sum for 
eq. ( p.20|) also has the same Lg dependence. The propagators in the 2- and 



3- direction have the same form. In the 4-direction Lg goes to Lf. 

Naturally the eigenvalues of the transfer matrix are independent of the 
temporal extension of the lattice, nevertheless the rate at which correlation 
functions decay can depend on Lj, and it is this decay rate which is used to 
measure ^r at finite Lt- The infinite volume limit turns out to be independent 
of Lt/Ls. 

The change in the solution of ( ^.141 ) will be proportional to the difference 
between the finite volume sum and the infinite volume integral when this 
difference is itself small. This leads to the expectation that for Lt > Lg 



3 



;Ur - yUR(oo) oc goilJ-R/Ls) 2 exp(-/iRLs) (3.23) 

and 



3 



ixx) - (xx)(oo) oc (^r/L,) 2 exp(-^RL,) = A . (3.24) 

These equations are to be understood as giving the form of the Lg depen- 
dence at fixed bare parameters for large Lg. The constant of proportionality 
depends on the values of the bare parameters. 

The above predictions have been tested by plotting {xx) or /iR against 
the right-hand side of (|3.23| ), (|3.24|) respectively. Data for a given lattice 



shape should fall on a straight line. In fig. ^ we show the Schwinger-Dyson 
results for the chiral condensate and see that the asymptotic formula ( p.24| ) 
holds for Lg > 6. A similar conclusion holds for the Monte Carlo results as 



fig. |T0| reveals. We gain added faith in this extrapolation to infinite volume 
from the fact that lattices with Lf = Lg and with Lf = 2Ls both extrapolate 
to the same value at infinite volumes. The extrapolated values for (xx) ^^^ 
/xr are given in tables |I], |^, |^ and ||. It is pleasing to note that the results 
from the 12"^ lattice already lie very close to the infinite volume results. 

The above formulae are only valid in a domain where any mass x Lg is 
large. Near the critical line 771^^ x Lg is very small and the finite size formulae 



of |22| should apply. Near the critical point all masses x Lg are very small 



and the exponential corrections turn into power law corrections. 
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4 Fermion-Antifermion Composite States 

We now turn to the determination of the energies of fermion-antiferniion 
composite states ("mesons"). 



4.1 Monte Carlo Method 

We calculated the correlation functions of bilinear operators of definite lattice 
symmetry which are well known from QCD spectroscopy (see, e. g., pi, 



We looked at four local operators and a one bond operator. For results 
on other one bond operators see |T^. The propagators for the local-local 
correlation functions have been calculated with 48 sources. Additionally we 
considered four wall operators |23|. The local operators O^^^ are given by 



Ol'^iit) =J2s,,k{x,t)x{x,t)x{x,t) 



(4.1) 



where the s factors and the corresponding continuum quantum number as- 
signments are given in the table |^. The notation for the continuum quantum 
numbers is the standard J^*-^. The subscript is the SU{4) flavour represen- 
tation (a for adjoint, s for singlet) under the assumption that the S'[/(4) 
flavour symmetry is restored. Our flfth operator is 



C'SW = -J2Vk{x,t)x{x,t) [x{x + k,t) + x{x-k,t) 



(4.2) 



For the wall operators we take 



l(w). 



Otkit) =( — ) J2 s^A^, tMx + e, t)xix, t) 



(4.3) 



i 


Si,kix,t) 


jPC 


name 


jPC 


name 


local operators 


1 


i-iy 


0.--*- 


vr' 


or 


a 


2 


/'_]^-\xi+a;2+X3+t 


o„-+ 


7r 


or 


— 


3 


C_na;fc+t 


K- 


P' 


ir 


a 


4 


/'_]_-\a;i+X2+X3-a;fe+t 


la" 


P 


ir 


h 


one-bond operator 




Vkix,t) 


1;^ 


to 


ir 


h' 


terms in eq. ( 


4:Ti|): 


E+ 


E- 



Table 7: Sign factors, corresponding continuum quantum number assign- 
ments and names of the particles that are used in the text. The last line 
gives the connection to the fit function, eq. (\4.11\ ). 
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The e sum extends over the vectors with all coordinates even. Correlation 
functions are 

Cf>it) = (OS(t)oS(0))-(Og(t))(o£(0)), / = 0,1, (4.4) 

Ct\t) = {Oi%t)O^J,\0))-{Ol%t)){Oij\0)). (4.5) 

They are independent of k. Whenever several values of k are possible we 
averaged over them. 

These correlation functions are the sum of a fermion line connected and 
a fermion line disconnected part (annihilation part). Pictorially this is rep- 
resented as 




O O - O O ^ (-) 



c 



where the brackets indicate averaging over the auxiliary field. Usually one 
measures only the fermion line connected part which is the first term in (^4.6|) . 
The second term, which we call the annihilation part, is computationally very 
demanding because the noise problem is more severe. However we have also 
attempted to measure the annihilation part with a stochastic estimator. The 
last term only contributes for particles with vacuum quantum numbers. 

Let us now turn to the discussion of our fit formulae. By means of a 
complete set of eigenstates of the transfer matrix one can write 

c{t) = {o{t)o'm 

= $:(z|a|j)(j|a'K)e-^^*-^»(^-*V*a^-*, (4.7) 

with (Tj = ±1. For a more detailed discussion of the O and the cxj see [P^] . 
If Lt was infinite one of the intermediate states in the sum would always be 
the vacuum, but on a finite lattice we have a genuine double sum. Terms 
will occur in the sum where Ei = Ej which contribute either a constant or 
an "oscillating constant" depending on whether a"j and aj have the same or 
opposite sign. (The reader may at first be surprised to hear that there are 
states with identical energy but opposite a. In the fermionic sector there 
are several discrete symmetries which change a but leave E unchanged. For 



example the action (|2.2|) is invariant under the transformation 

X{x) -^ C^{x)x{x + A), x{x) -^ C^{x)x{x + /i), (4.8) 

where 

This transformation changes a if /i is spacelike. Such pairs occur in the 
fermionic sector because of the well-known phenomenon of doubling. The 
propagator for staggered fermions has sixteen states instead of the four states 
of the continuum propagator. Eight of these states come from poles with p^ 
near 0, and so correspond to ordinary exponentials in position space, the 
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Figure 11: Meson propagators Ci{t) from Monte Carlo calculations on an 
8^ X 16 lattice at mo = 0.02 and go = 0.21 (l.h.s.), go = 0.32 (r.h.s.). Shown 
are the local-local propagators C^ (t) (circles) and the wall-local propaga- 
tors Cl (t) (diamonds). Open (solid) symbols represent positive (negative) 
values. The lines connect values from simultaneous fits using eq. (\4.1]\) . 



24 



,,(0) 


90 


12^ 


8^ X 16 


mo = 0.01 


mo = 0.02 


mo = 0.02 mo = 0.04 


mo = 0.09 


0.21 


0.095(4) 




0.191(1) 


0.326(1) 


0.582(1) 


0.23 


0.129(4) 




0.238(1) 


0.377(1) 




0.25 


0.179(4) 


0.298(2) 


0.280(1) 


0.420(1) 




0.26 






0.297(1) 


0.434(1) 




0.27 


0.238(5) 


0.344(3) 


0.310(1) 


0.442(1) 


0.653(1) 


0.28 


0.250(3) 




0.317(1) 


0.445(1) 




0.29 






0.319(1) 


0.444(1) 


0.644(1) 


0.30 


0.244(2) 


0.334(2) 


0.317(1) 


0.442(1) 


0.638(1) 


0.32 


0.232(1) 




0.311(1) 


0.431(1) 


0.624(1) 



,,(1) 


9o 


12^ 


8^ X 16 


mo = 0.01 


mo = 0.02 


mo = 0.02 


mo = 0.04 


mo = 0.09 


0.21 


0.772(5) 




0.880(7) 


0.923(6) 


1.024(5) 


0.23 


0.678(7) 




0.828(10) 


0.885(8) 




0.25 


0.594(11) 


0.646(8) 


0.807(10) 


0.889(16) 




0.26 






0.816(14) 


0.885(11) 




0.27 


0.570(24) 


0.669(13) 


0.831(21) 


0.947(20) 


1.124(15) 


0.28 


0.635(20) 




0.866(21) 


0.974(15) 




0.29 






0.902(27) 


1.014(11) 


1.232(15) 


0.30 


0.800(27) 


0.885(32) 


0.974(33) 


1.064(15) 


1.273(21) 


0.32 


0.953(40) 




1.096(35) 


1.201(33) 


1.375(22) 



Table 8: Monte Carlo results for the energies of the ground state of the vr 
particle E^^ and the first excited state E^^ . 



other eight have p^ near to vr and therefore give rise to oscillating expo- 
nentials. The above symmetries ensure that the decay rate of the normal 
exponentials and the oscillating exponentials remains the same, even in the 
presence of interactions. Therefore we know that every single fermion state 
has a partner state with the same energy but opposite a, so there will indeed 
be contributions of the "oscillating constant" type in the meson propagator. 
We have already come across a pair of states with the same energy but op- 
posite a in the fermion propagator eq. ( p.l5| ), fig. |^ The reader can easily 
check that at small t the propagator G{t) is dominated by a normal expo- 
nential and at large t by an oscillating, but that in both cases the energy is 
the same.) Consider two processes contributing to the correlation functions: 

(^) "ClI^ • The intermediate state is a fermion-antifermion pair, 

(b) -^'''''^N-i- : A single fermion propagates around the lattice. 

It is process (b) which gives rise to the constants. This can already be 
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checked for the free case on a finite lattice. The check is easiest for the wall 
propagator C^^™ which in the free theory is simply the square of the fermion 
propagator G{t) (see eq. ( p.l5| )) 

Ci^^{t) = GlJt) oc e^2^'^* + e~2^^(^*~*) - 2(-l)*e-^«^* (4.10) 

which already displays an oscillating constant contribution. 
From eq. ( [4.7| ) C{t) thus has the form 

C{t) = $:yl+(e-^"* + e-^"(^*-*)) 

n 

+ E^n (e""""* + e-^"(^*-*)) (-1)* 

n 
+ 5+ 

+ B-{-iy. (4.11) 

The E~^ terms describe s-wave states and the E~ terms describe the parity 
partners of these states which are p-waves. The quantum numbers are given 
in table |^. In the following the energies are called -E^ame where "name" is 
the particle name from table ^ and i = or 1 for the energies of the ground 
state or first excited state. 

In fitting we had to restrict ourselves to the lightest states. We used two 
s-wave states for the vr and one s- and one p-wave state for the tt'/ct. For all 
vector channels {p' /a, p/b, uj/b') we used two s-wave states and one p-wave 
state. The choices of the terms in the fit function were necessary to obtain 
good fits. In the case of the vr correlation function it is further motivated 
by the fact that a O"*"" state cannot be realised as a fermion- ant if ermion 
state. In the other cases we were guided by the results from the Schwinger- 
Dyson equations. It was found that B^ was always negligible and that B~ 
was necessary to describe propagators in the symmetric phase. In order to 
extract a maximum of information from our correlation functions we found 
it convenient to make a simultaneous fit to both the local-local and the wall- 
local correlation functions. 

From the vr ground state energy E^^ and amplitude A^^ we calculated 



the pion decay constant /^^ using p5|]: 



mo 



4 



l(0) 



sinh E: 



(0) 



f^ = -^\^'^' ■\..^rn ,.. ■ (4.12) 



sinh'(Ei°V2) 



The fit results are given in tables P- 13. As we have seen for the chiral conden- 



sate the finite size effects are quite large. In fig. |TI| we show propagators and 
fits of the vr'/cr, vr and p/b (top to bottom) for one gQ-valne in the symmetric 
and one Qo-vaAue in the broken phase. 
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Table 9: Monte Carlo results for the pion decay constant /^ 
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1.032(8) 


0.28 


0.501(25) 




0.675(26) 


0.830(15) 




0.29 






0.747(33) 


0.907(11) 


1.130(12) 


0.30 


0.683(53) 


0.783(25) 


0.796(44) 


0.981(28) 


1.223(16) 


0.32 


0.892(110) 




0.934(72) 


1.131(45) 


1.303(24) 



,,(0) 


90 


12^ 


8^ X 16 


mo = 0.01 


mo = 0.02 


mo = 0.02 


mo = 0.04 


mo = 0.09 


0.21 


0.751(8) 




0.877(6) 


0.939(4) 


1.105(3) 


0.23 


0.665(9) 




0.817(6) 


0.922(4) 




0.25 


0.595(11) 


0.680(6) 


0.811(5) 


0.936(9) 




0.26 






0.816(7) 


0.940(5) 




0.27 


0.598(17) 


0.705(11) 


0.824(11) 


0.974(12) 


1.206(10) 


0.28 


0.656(15) 




0.855(11) 


0.996(8) 




0.29 






0.854(14) 


1.026(6) 


1.271(9) 


0.30 


0.771(25) 


0.873(15) 


0.880(17) 


1.049(23) 


1.300(19) 


0.32 


0.898(47) 




0.926(29) 


1.152(21) 


1.364(17) 



Table 10: Monte Carlo results for the energies of ground state of the ir' 
particle E^, and the ground state energy of the a particle E^\ 
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^(0) 


90 


12^ 


8^ X 16 


mo = 0.01 


mo = 0.02 


mo = 0.02 


mo = 0.04 


mo = 0.09 


0.21 


0.101(2) 




0.253(1) 


0.407(1) 


0.678(1) 


0.23 


0.144(2) 




0.344(2) 


0.499(2) 




0.25 


0.227(3) 


0.371(2) 


0.460(3) 


0.620(2) 




0.26 






0.526(4) 


0.677(3) 




0.27 


0.374(5) 


0.521(4) 


0.595(5) 


0.749(4) 


0.992(4) 


0.28 


0.485(5) 




0.662(7) 


0.817(6) 




0.29 






0.741(9) 


0.883(4) 


1.111(5) 


0.30 


0.676(11) 


0.776(7) 


0.829(13) 


0.935(11) 


1.187(11) 


0.32 


0.890(29) 




0.969(43) 


1.076(17) 


1.290(19) 



i,(i) 


90 


12^ 


83 X 16 


mo = 0.01 


mo = 0.02 


mo = 0.02 


mo = 0.04 


mo = 0.09 


0.21 


1.36(1) 




1.40(2) 


1.46(2) 


1.54(1) 


0.23 


1.34(1) 




1.42(2) 


1.48(2) 




0.25 


1.30(2) 


1.35(1) 


1.47(2) 


1.47(4) 




0.26 






1.38(3) 


1.43(2) 




0.27 


1.34(4) 


1.34(3) 


1.39(4) 


1.61(7) 


1.56(3) 


0.28 


1.52(7) 




1.52(7) 


1.44(3) 




0.29 






1.42(6) 


1.43(2) 


1.57(3) 


0.30 


1.57(11) 


1.46(5) 


1.52(9) 


1.51(10) 


1.58(6) 


0.32 


1.32(10) 




1.21(9) 


1.60(10) 


1.53(5) 



Er 


90 


12^ 


8^x 16 


mo = 0.01 


mo = 0.02 


mo = 0.02 


mo = 0.04 


mo = 0.09 


0.21 


1.41(3) 




1.49(5) 


1.65(5) 


1.73(3) 


0.23 


1.42(4) 




1.48(6) 


1.69(6) 




0.25 


1.35(5) 


1.46(4) 


1.63(7) 


1.64(11) 




0.26 






1.52(8) 


1.63(6) 




0.27 


1.53(15) 


1.40(10) 


1.51(13) 


1.97(20) 


1.77(11) 


0.28 


1.98(24) 




1.97(22) 


1.66(9) 




0.29 






1.55(15) 


1.62(6) 


1.75(8) 


0.30 


1.90(31) 


1.60(13) 


1.72(21) 


1.99(34) 


1.77(14) 


0.32 


1.23(19) 




1.45(21) 


2.11(29) 


1.82(11) 



Table 1 1 : Monte Carlo results for the energies of the ground state of the p 
particle EJ^\ its first exited state EJ^^^ and the ground state of the b particle 



E^ 



(0) 
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p(0) 


90 


12^ 


8^x 16 


mo = 0.01 


mo = 0.02 


mo = 0.02 


mo = 0.04 


mo = 0.09 


0.21 


0.100(2) 




0.263(2) 


0.418(2) 


0.692(2) 


0.23 


0.140(2) 




0.352(3) 


0.509(3) 




0.25 


0.224(3) 


0.374(2) 


0.470(4) 


0.633(3) 




0.26 






0.534(5) 


0.688(4) 




0.27 


0.380(6) 


0.519(5) 


0.608(7) 


0.770(5) 


1.016(5) 


0.28 


0.490(7) 




0.680(10) 


0.830(7) 




0.29 






0.767(12) 


0.904(5) 


1.141(7) 


0.30 


0.685(14) 


0.784(10) 


0.837(16) 


0.943(12) 


1.214(13) 


0.32 


0.840(25) 




0.930(31) 


1.099(17) 


1.282(11) 





90 


124 


8^x 16 


mo = 0.01 


mo = 0.02 


mo = 0.02 


mo = 0.04 


mo = 0.09 


0.21 


1.40(3) 




1.46(5) 


1.61(4) 


1.68(2) 


0.23 


1.41(4) 




1.48(6) 


1.61(5) 




0.25 


1.40(5) 


1.44(4) 


1.64(7) 


1.60(9) 




0.26 






1.66(10) 


1.63(5) 




0.27 


1.48(8) 


1.64(10) 


1.51(13) 


1.76(13) 


1.80(8) 


0.28 


1.50(12) 




1.51(12) 


1.61(7) 




0.29 






2.08(33) 


1.93(9) 


1.90(7) 


0.30 


1.78(26) 


1.58(11) 


2.04(34) 


1.70(20) 


1.68(11) 


0.32 


2.07(57) 




1.70(27) 


2.12(31) 


2.23(18) 



E^^^ 


9o 


12^ 


8^ X 16 


mo = 0.01 


mo = 0.02 


mo = 0.02 


mo = 0.04 


mo = 0.09 


0.21 


1.36(1) 




1.43(2) 


1.48(1) 


1.57(1) 


0.23 


1.34(1) 




1.42(2) 


1.50(2) 




0.25 


1.32(2) 


1.36(1) 


1.53(2) 


1.50(3) 




0.26 






1.49(3) 


1.49(2) 




0.27 


1.38(3) 


1.46(4) 


1.46(4) 


1.61(5) 


1.71(3) 


0.28 


1.40(4) 




1.46(4) 


1.54(3) 




0.29 






1.75(11) 


1.65(3) 


1.75(3) 


0.30 


1.51(8) 


1.48(4) 


1.66(10) 


1.61(8) 


1.70(4) 


0.32 


1.79(22) 




1.49(8) 


1.79(10) 


1.88(6) 



Table 12: Monte Carlo results for the energies of the ground state of the p' 
particle E , , its first exited state E , and the ground state energy of the a 
particle EJ^^ . 
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E^S^ 


9o 


12^ 


S^x 16 


mo = 0.01 


mo = 0.02 


mo = 0.02 


mo = 0.04 


mo = 0.09 


0.21 


0.22(4) 




0.27(1) 


0.43(1) 


0.71(1) 


0.23 


0.30(4) 




0.37(2) 


0.53(1) 


0.80(2) 


0.25 


0.46(5) 


0.47(4) 


0.50(4) 


0.64(3) 


0.94(4) 


0.26 






0.62(4) 


0.71(4) 


1.00(6) 


0.27 


0.56(13) 


0.40(15) 


0.82(8) 


0.78(5) 


1.11(5) 


0.28 


0.54(12) 




0.58(9) 


0.80(8) 


1.01(14) 


0.29 






0.80(12) 


1.04(7) 


1.32(8) 


0.30 


0.99(10) 


1.03(4) 


1.29(9) 


0.91(21) 


1.29(15) 


0.32 


1.22(11) 




0.97(20) 


0.97(15) 


1.18(25) 



(jj 


90 


12^ 


8^ X 16 


mo = 0.01 


mo = 0.02 


mo = 0.02 


mo = 0.04 


mo = 0.09 


0.21 


1.68(3) 




1.75(4) 


1.76(3) 


1.78(2) 


0.23 


1.74(6) 




1.68(7) 


1.75(5) 


1.81(3) 


0.25 


1.97(18) 


1.72(6) 


1.79(18) 


1.78(6) 


1.90(6) 


0.26 






1.86(14) 


1.86(11) 


1.90(7) 


0.27 


1.82(34) 


1.57(10) 


2.23(51) 


1.85(11) 


2.05(11) 


0.28 


1.68(21) 




1.64(10) 


1.80(9) 


1.88(9) 


0.29 






1.80(22) 


2.22(27) 


2.28(27) 


0.30 








1.74(12) 


2.12(25) 


0.32 








2.13(31) 


1.98(15) 



Eb' 


9o 


12^ 


8^ X 16 


mo = 0.01 


mo = 0.02 


mo = 0.02 


mo = 0.04 


mo = 0.09 


0.21 


1.55(4) 




1.67(7) 


1.75(5) 


1.71(3) 


0.23 


1.67(7) 




1.42(18) 


1.64(13) 


1.67(8) 


0.25 


1.90(18) 


1.64(7) 


1.83(48) 


1.70(8) 


1.79(11) 


0.26 






1.57(23) 


1.98(30) 


1.91(16) 


0.27 


1.69(34) 


1.56(14) 


2.15(46) 


1.82(13) 


1.93(10) 


0.28 


1.49(21) 




1.27(26) 


1.68(23) 


1.70(16) 


0.29 






1.69(24) 


2.15(22) 


2.15(19) 


0.30 








1.61(14) 


2.04(17) 


0.32 








2.33(46) 


2.02(19) 



Table 13: Monte Carlo results for the energies of ground state of the u 
particle E^^ , its first exited state E^"* and the ground state of the b' particle 



E 



(0) 



30 



C2(t) 




C2(t) 



10 - 




Figure 12: Influence of tie anniliilation part on the n propagator in tie 
Monte Carlo calculation. Parameters and symbols are the same as in fig. [7I| 
except that the diamonds now represent the complete propagator, eq. ( |4.6| ). 





Figure 13: Influence of the annihilation part on the vr'/cr propagator in the 
Monte Carlo calculation. Parameters and symbols are the same as in fig. [7I| 
except that the diamonds now represent the complete propagator, eq. ij^ 
See text for comments on (b). 
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For certain values of Qq and rriQ we have also measured the annihilation 
part of the vr and a correlation functions, where 16 sets of gaussian random 
numbers were used in the stochastic estimator. As fig. |1^ shows the influence 
in the tt correlation function is small. Since the a has the quantum numbers of 
the vacuum we have to subtract a constant from the a correlation function 
before comparing it with the result without the annihilation part. As the 
constant is rather large (it is of 0(L^(xx)^)) i^i the middle of the correlation 
function in the broken phase we only see noise. A typical result is given in 
fig. O. Within the errors the energies are unchanged and thus we can use 



the numbers in table 10 with some confidence. The discussion of our results 



will be postponed until section ^ 



4.2 Schwinger-Dyson Equations 

The fact that the Schwinger-Dyson equations have given such good results for 
the fermion propagator and chiral condensate encourages us to try a similar 
approach to meson propagators. In fig. [l^ we show the Schwinger-Dyson 
equations for the fermion-antifermion composite propagator. 

To get a manageable set of equations we have replaced the full kernel with 
the bare 4-fermi interaction. This gives us equations which are at the same 
level of accuracy as the order go Schwinger-Dyson equations for the fermion 
propagator which, as we have already seen, yield a good approximation to 
the Monte Carlo results. 





M 



+ 




M K P 



Figure 14: Pictorial representation of the Schwinger-Dyson equation for the 
meson propagators. 



32 



For the ferinion line connected part of the local correlation functions we 
derive the following results in appendix 0. Defining 






(0), 



we find that the Fourier transforms C\ (cu) are given by 



(4.13) 



where 



Cf\u) 



ai{uj) = — 



'e 

v^ 



ai{uj) 



1 + ai{uj)2go{+3 - 


- COSCj) 


a2{uj) 




1 + a2{uj)2go{-3 - 


- COS a;) 


a-iiuj) 




l + a3H2^o(+l- 


- coscj) 


04(0;) 





1 + a4{io)2gQ{—l — coscj) 



Bf + F^ sin ^4 sin(fc4 — tu) 



[A2 + Ft^ sin^ k^ 



A'^ + Ft^ sm\h - u] 



(4.14) 



(4.15) 



with 



and 



A' 


= 


A^' + F^s' , 


Bl 


= 


N' - F^s' , 


Bl 


= 


N^ + F>2 , 


Bl 


= 


N^ - -Fh^ 
3 ' 


Bl 


= 


N' + \Fy 



s^ = ^ sin^ ki . 



i=l 



(4.16) 



(4.17) 



The iV, Fs = Fi = F2 = F3 and Ft = F4 have been defined in (^J§. The 
formulae for the fermion line disconnected part are also given in appendix ^ 
The principal interest in the fermion- ant if ermion channels lies in the en- 
ergy spectrum, i. e., the eigenvalues of the transfer matrix. These energies 
can best be found on a long lattice which is why we used 8^ x 16 lattices 
for most of our simulations. In the Schwinger-Dyson approach we can easily 
extend Lt to infinity, which allows a clean extraction of the full spectrum. 
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When Lt = oo the ^4 sum becomes an integral, which can be evaluated 
by contour integration. The general "bubble" Oj in (|4.15| ) has the form 



«^M = ^T.AA^B,u), (4.18) 

* k 

IiA,B,u) = r^, B^ + F,^ sin hsu^ih-uj) 

Jo 271 ^A^ + Fism^k4\[A^ + Ft^sm^{h-uf 

where A and B are functions of k only. Evaluating this integral leads to the 
result 

1 / A^ + B^ A^-B^ \ 

^ ' ' ""^ ~ aJa^ + F^ U^' + Ft - F^ cos a; 2A^ + F^ + F^ cosij) ' 

(4.20) 
From the form of X it is easy to see that the Fourier transform of the 
composite propagator is the ratio of two real polynomials in coscj. There- 
fore the only singularities that can occur when Lg is finite are poles. A 
pole in the momentum-space propagator at u = iE corresponds to a term 
oc exp{—Et) in the real-space propagator, and a pole aX u = tt + iE to a 
term oc (— l)*exp(— ii^t), so we can find the complete energy spectrum by lo- 
cating all the poles in the momentum-space propagator. The order of these 
polynomials is proportional to the number of distinct terms in the sum over 
spatial momenta, so the number of energy levels should grow as L^. 

Cursory inspection of the meson levels in tables |8|-p!3| shows that many 
of the energy levels lie close to or above 2/iR. In order to really understand 
what is happening above the fermion-antifermion threshold we need to take 
the infinite volume limit and see what happens as the excited states form a 



true continuum. We have already taken the Lt ^ 00 limit in eq. ( [4.20|) . We 
now let Ls —^ 00. The sum over k in eq. ( [4.15| ) then turns into an integral, 
the analytic structure of which is worked out in appendix ^. As Lg increases 
the poles above the threshold become denser and in the infinite volume limit 
become a cut. 



In fig. Il5] we sketch the singularity structure of C in the cj-plane and 



also give the integration contour needed to evaluate C{t). C{uj) has cuts 
running from cos cj = 1 + 2 sinh /ir to cos u; = 7 + 2 sinh /xr and from 
cosu; = — 1 — 2sinh^/iR to cosu; = —7 — 2sinh^/iR. These cuts represent 
the continuum of fermion-antifermion states running from E = 2/iR to ii^ = 
cosh~^(7 + 2smh.^ fx^), the highest energy available to two fermions. The 
singularities in the phase-space function (f){r) (see appendix B) at r = ±1 
lead to additional branch points within the cuts at coscj = ±(3 + 2 sinh^ /xr) 
and cos cj = ± (5 + 2 sinh^ /ir) . 



The discontinuity across the cuts of C gives the spectral function p [^ 
As well as these cuts C can also exhibit poles, which correspond to genuine 
bound states (stable mesons). 
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CJ 





7T 



1 



1 






1 



Figure 15: Singularity structure of C, 



(0), 



UJ 



in the uj-plane. The sohd hnes 



are cuts, a * is a pole and a • is a branch point. 



4.3 Finite Size Effects 

We shall now study finite size effects for bound state masses using methods 



similar to those used in section |3.3| for the chiral condensate and the fermion 
mass. To do this we will consider finite size corrections to the "bubble" 
sums given in eq. ( |A.7D . The main characteristic of these bubbles is that 



they involve two propagators with denominators U^ and U~ . We take a 
slightly more general case where the physical masses (m+ and m_) in the 
two propagators can be different. 

To find the change in the meson propagator we must investigate the 
effect of replacing the sum over loop momenta by an integral. Since we are 
interested in the meson mass we need to consider the case where the meson 
is on-shell, i. e., the 4-momentum p = (0, 0, 0, iM) where M is the meson 
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mass. We now use the same Fourier transform argument as in section p.3| to 
change a sum in momentum space to a sum in coordinate space. 
A typical loop integral in eq. ( |A.7| ) has the form 

d^k n{p,k) 

(27r)4((p + A;)2 + m2.)(P + m^) ' ^ ' ^ 

(The exact form of the numerator n{p, k) is unimportant for asymptotic 
behaviour.) As before we need the large Lg limit of 



d^k 



exp{ikiLs 



n{p, k) 



(4.22) 



(27r)4 --^v— «/ ((p + ky + ^2_) (p ^ ^^) ■ 

Integrating exactly over ki and then with the saddle point approximation 
over the other three components of k gives a result 



oc — exp(-^L,) = A, 

L/s, 



(4.23) 



where 



e 



-M4 + 2M^ml + 2M^ml + 2mlm^ - mi - mt 



2M 



(4.24) 



The finite-size shift in the meson mass ought to be proportional to A. 




Figure 16: Dispersion relation of the ii bound state obtained from Schwin- 
ger-Dyson calculations. Results are from a 52^ x oo lattice with mo = 0.01 
and go = 0.2476, leading to fiR = 0.20 and nij, = 0.301. 
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We can consider various special cases of eq. (^4.24 ) . In the non-relativistic 



limit where the binding energy B = m^ + m_ — M is small compared with 

M 



C ~^ V 2"^reduced-B 



^^^ ^ {m+ + m_- M), (4.25) 



\ \"^+ + Tn 



recovering a result that can be found by solving the Schrodinger equation in 
a finite volume. 

Other useful limits are the case m+ = m^ which gives ^ = {Am^ — 
M^)^/^/2, which should apply to the pion when it is a true bound state, and 
M = m+ = ru- which gives ^ = ^M, relevant for finding the pion mass 
shift due to tt ^ 27r ^ vr [^. A further application of eq. ( [4.23 ) is to the 



process / — >■ Z + tt — i> /, which will give a finite size shift to the fermion mass 



IXr additional to that considered in section 3.3. This extra contribution will 
be of the form ( |4.23| ) with ^ = ^^ v 4/u|^ — m^ and can become the leading 



term if the pion is light enough. 

With the Schwinger-Dyson equations it is easy to study the restoration 
of Lorentz invariance. As an example in fig. |l^ we look at the dispersion 
relation of the vr bound state. We have plotted E"^ against p ^ for all possible 
p values which occur on a lattice with Lg = 52. When p^ can be realised 
in inequivalent ways we have plotted all the possible energies. We see that 
the Einstein relationship E"^ = rn^ + p ^ holds well up to p"^ k, Q.2 and that 
spherical symmetry [E'^ independent of the momentum direction) holds up 
to p^ ^ 0.5. 

Although Lorentz invariance is restored we observe at most partial flavour 
symmetry restoration: The correlation functions for the Goldstone vr and the 
non-Goldstone vr look rather different. 
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5 Meson Propagators and Spectroscopy 



We now compare the correlation functions found by solving the order qq 
Schwinger-Dyson equations with those from Monte Carlo. We show some 
Monte Carlo propagators together with the Schwinger-Dyson correlation 
functions calculated with the same bare mass and the coupling chosen such 
that the renormalised fermion masses agree. In view of the fact that the order 
go Schwinger-Dyson equations do not reproduce exactly the critical coupling 
found in the Monte Carlo calculations it is not surprising that a shift in the 
coupling is necessary to achieve agreement in the propagators. Indeed the 
shift is about the same as the difference between the critical couplings. Look- 



ing at figures p!7H19| we find a good agreement between the Schwinger-Dyson 
and Monte Carlo data. The agreement is of course best at small Qq but even 
at large go is still satisfactory. 



C2(t) 




Figure 17: Comparison of Schwinger-Dyson and Monte Carlo data for the vr 
propagators C2 (t) (sohd hnes, circles) and C2 (t) (dashed lines, diamonds). 
The results come from calculations on an 8^ x 16 lattice at mo = 0.02 and 
go = 0.27 (Monte Carlo). For the Schwinger-Dyson calculation we used 
go = 0.2491 leading to the same fermion mass as the Monte Carlo calculation. 
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Ci(t) 




Figure 18: Comparison of Schwinger- Dyson and Monte Carlo data for the 
7['/(j propagators. Symbols and parameters are the same as in fig. |I^. Solid 
symbols represent negative values. 



Like the Monte Carlo measurements the Schwinger-Dyson equations give 
a rather small annihilation contribution to the correlation functions. The 
largest effect is seen for the a particle deep in the broken phase. An example 



is shown in fig. ^ In distinction to the Monte Carlo approach there are no 
noise problems. Let us note that the order go Schwinger-Dyson equations 
satisfy the Ward identity EB, m 



d{xx) 

dlTLn 



tr^(o). 



E(-i)*cr(i) 



(5.1) 



if one uses the full correlation functions (but not if the annihilation contri- 
bution is neglected). 

From now on we just look at the fermion line connected part of the 
propagators. In figures pT|-p^ we compare energy levels from the Monte 
Carlo with those from the Schwinger-Dyson calculations. In these figures we 
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C4(t) 




Figure 19: Comparison of Schwinger- Dyson and Monte Carlo data for the 
p/b propagators. Symbols and parameters are the same as in fig. |I^. Solid 
symbols represent negative values. 



plot some measured values from tables |8|-pT] against the renormalised fermion 
mass. The s-wave states are represented by solid lines (Schwinger-Dyson) and 
solid circles (Monte Carlo), the p-wave states by dashed lines (Schwinger- 
Dyson) and open circles (Monte Carlo). The dotted lines represent 2/iR, so 
bound states lie below this line. The states that lie above the threshold would 
form the continuum on the infinite lattice possibly including resonances. One 
must be careful about states that lie slightly below the threshold because they 
might appear to be bound due to the finite volume. 

Let us discuss the pictures in more detail. The vr picture (fig. ^l]) shows 
that the first two levels are in very good agreement. It is thus reasonable 
to expect that the higher levels are correctly represented by the Schwinger- 
Dyson formalism. The large gap between the second and third level explains 
the success of the two level fit formula. Obviously we find a bound state 
with the quantum numbers of the vr in the broken phase. This state is the 
pseudo-Goldstone boson associated with chiral symmetry breaking, and its 
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Figure 20: Influence of tlie anniliilation part on tlie n' jo propagator in tlie 
Scliwinger-Dyson calculation. The solid line is CJ (t) calculated without the 
annihilation part, the dashed line represents the result including it. The 
parameters are the same as in fig. \l^(b), except that go = 0.2975 leading to 
the same fermion mass as the Monte Carlo data at Qq = 0.32. 



mass tends to zero as mo tends to zero as one expects. The picture shows 
avoided level crossing suggesting the possible existence of a resonance [|3 ^^ 
the symmetric phase. 



In the vr'/o" picture (fig. |22|) we see two states with opposite parity (the 
a is the top curve). For small /ir the ground state energies are reproduced 
by the Schwinger- Dyson equations. There is somewhat less good agreement 
for large values of /ir. It is not clear whether we have a slightly bound state 
for the a. The ground state energy of the n' is always close to 2/iR. 

Finally in fig. |2^ we show levels in the p/b channel (the p is the solid 
curve). The agreement for the p ground state is good, for the excited states 
the fit cannot really resolve the levels lying relatively close to each other. 
The ground state of the p does not appear to be bound. The p' /a and u/b' 
pictures look similar. 
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Figure 21: Energy levels in the it channel on 8^ x 16 (Monte Carlo) and 
8^' X oo lattices (Schwinger-Dyson) at itlq = 0.02. The dotted line represents 
the threshold Ej, = 2/iR. 
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Figure 22: Energy levels in the it' /a channel on 8^ x 16 (Monte Carlo) and 
8^ X oo lattices (Schwinger-Dyson) at ttiq = 0.02. The dotted line represents 
the threshold E^^r/^ = 2/iR. The solid line and the black symbols represent 
E^i while the dashed lines and open symbols represent E^. 
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Figure 23: Energy levels in the p/b channel on 8^ x 16 (Monte Carlo) and 
8^ X cx) lattices (Schwinger-Dyson) at mo = 0.02. The dotted line represents 
the threshold Ep/b = 2/iR. The solid line and the black symbols represent Ep 
while the dashed lines and open symbols represent Ej,. 



We now look at what is happening in large volumes. In fig. ^ we have 
plotted the energy levels of the vr on a 52^ x oo lattice (solid lines). Also 
shown are the ground state energies for smaller lattices together with Monte 
Carlo data. For large /xr we have a bound state which for small values of /ir 
turns into a resonance, whose energy was calculated from the infinite volume 
spectral function (see below). In the finite volume the resonance reveals itself 
by avoided level crossing. Finite size effects for the ground state are largest 
close to the threshold. 

For comparison we also display energy levels of the a and p on a large 
lattice along with the resonance energies (figures ^ and ^6]). In both cases 
there are bound states for large /zr. We would like to point out that the p 
resonance and bound state energies stay large even as mo — ^ (see below). 
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Figure 24: Energy levels of the vr particle on a 52^ x oo lattice (mo = 0.02). 
Also shown are the ground state energies on 20^ x oo, 12^ x cxd and 8'^ x cxo 
lattices together with Monte Carlo results from 12'* and 8^ x 16 lattices. The 
solid line indicates the resonance energy at infinite volume. 



As the volume increases the density of levels above threshold increases 
further until for infinite volume we have a continuum and can now define 
spectral functions p^, / G {0, 1, w}: 



dE 

IT 



C(0(t) = / !^e~^l*lp?(^) + (-1)* / ^e-^'IV^'(^) • (5.2) 



^t / C^-S^-Bltl (0, 



TT 





They are calculated from the discontinuities across the cuts and poles of the 



amplitude C in fig. |T5: 

p\ = lm&{iE + e) 



P- 



\mCHiE + Tx + t) , e 



-0. 



(5.3) 
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Figure 25: Energy levels of the a particle on a 20'^ x oo lattice (mo = 0.02J. 
The solid line indicates the resonance energy at infinite volume. The dotted 
line is the threshold. 

Let us look at the vr spectral function in more detail. In this case we 
need only one spectral function /9+, because p_ is zero. Above 2/iR there 
is a continuum of states (fig. P7[b). In the symmetric phase fig. P7| a shows 
a resonance. We can tell that this 'bump' is really a resonance because we 
have analytically continued the amplitude to the second Riemann sheet and 
found a pole. Locating the pole gives us both the real and imaginary parts 
of the pion mass and so tells us the rate T^^ at which the tt resonance decays 
to a fermion-antifermion pair. In the broken phase the spectral function of 
the IT contains a (5-function contribution. This shows that the resonance has 
turned into a bound state. 



45 



E 



1.5 



0.5 




0.2 



0.4 



0.6 



0.8 



Mr 



Figure 26: Energy levels of the p particle on a 20'^ x oo lattice (mo = 0.02J. 
The solid line indicates the resonance energy at infinite volume. 



In fig. ^ we show the vr resonance mass r/iTr and width F^r in the symmetric 
phase. Note that the ratio F^/m^r decreases towards zero as we approach the 
critical point, suggesting that the physical strength of the meson- fermion 
coupling decreases as we go to the continuum limit. This is what we would 
expect in a 'trivial' theory. 

In fig. ^ we show the spectral function of the a. There is a strong and 
narrow resonance just above the threshold. This corresponds to the scalar 
state at 2yUR that Nambu and Jona-Lasinio found M. 
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Figure 27: Spectral function of the n particle: (a) in the symmetric phase 
(m,Q = 0.04, go = 0.2056 leading to jir = 0.20J, (b) near gc (rriQ = 0.01, 
go = 0.2476 leading to fiR = 0.20;. 
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Figure 28: Resonance mass (solid line) and width (dotted line) of the % 
particle at niQ = 0. 




Figure 29: Spectral function of the a (dashed line) and vr' (solid line) particles 
near g^ (mo = 0.01, qq = 0.2476 leading to fiR = 0.20J. (Note that the two 
functions have been shown at very different scales.) 
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Figure 30: Spectral function of the p (solid line) and b (dashed line) particles 
in the symmetric phase (mo = 0.04, go = 0.2056, leading to pn = 0.20J. 
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Figure 31: Resonance mass (solid line) and width (dotted line) of the p 
particle at mo = 0.02. 
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Next we look at the picture of the p and h spectral functions (fig. ^0]). We 
see that the p is a resonance. Although this would initially appear interesting 
we find that this resonance never becomes light compared with the inverse 
lattice spacing. Even at the critical point it has a mass of about 1.5 (see 
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fig. ^). In the b channel there is only a continuum. Note that the s- and 
p-wave show a different threshold behaviour (p oc [-E — 2^^]^/^ for the s-wave, 
p (X [E — 2/iR]^/^ for the p-wave). 

As well as the resonance we see two other curious features, sharp discon- 
tinuities in the slope of the spectral function at energies of approximately 
1.8 and 2.3. In general they appear at ii^ = cosh~^(3 + 2sinh^pR) and 
cosh~^(5 + 2sinh^/iR). These are examples of van Hove singularities [P^]. 
Such singularities are familiar from solid state physics and result from saddle 
points in the fermions' energy momentum relation. At the corresponding en- 
ergies the density of states is singular leading to singularities in the spectral 
functions (see appendix 0). These singularities are also present (though often 
less noticeable) in the other channels. The van Hove singularities introduce 
a complication when analytically continuing onto the non-physical Riemann 
sheet (see fig. ^). Because there are branch points we reach different sheets 
depending on where we cross the axis in relation to these singularities. The 
sheet which is relevant for the continuum limit of the theory is the one reached 
by crossing the axis between the threshold and the first van Hove singular- 
ity. The other sheets only exist because of the saddle points in the lattice 
fermion dispersion relation at k = (|, 0, 0), (|, |, 0) etc., and so should have 
no relevance to the continuum. All the resonances we have seen are on this 
'low energy' sheet. 

The p' /a channel is similar to the p/b. We see similar resonances (one of 
each parity) which in the continuum limit go to the cutoff. For the u/b' the 
Monte Carlo results indicate that these particles are also heavy. 
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6 Renormalisation Group Flow 

In this section we study the renormahsation group flow in our model. We con- 
sider dimensionless ratios of physical quantities and attempt to find lines of 
constant physics, i. e., lines where all those ratios are constant, independent 
of the cutoff. Only in regions where such lines exist is the theory renormalis- 
able 1 11, 29, 3^. Since we are interested in physical quantities we should only 
look at the ratios involving the particles whose correlation lengths diverge at 
the critical point, namely the fermion, vr and a. Other states (such as the p) 
with energies near 1/a are irrelevant in this context. 

In fig. ^ we show Ef^/jiYi and /xr/Ztt for the Monte Carlo results on 
an 8'^ X 16 lattice. There is no region where the lines are parallel to each 
other, thus there are no lines of constant physics. There are however two 
caveats to this picture. The first is that we have large finite size effects and 
secondly that in the symmetric phase the tt is a resonance, and so Ef^ no 
longer corresponds to the pion mass (see fig. ^4]). The upper right hand 
corner suffers least from these problems, and here the Monte Carlo results 
can be taken at face value. These difficulties can be circumvented by use of 
Schwinger-Dyson equations. In fig. ^ we show the Schwinger-Dyson results 
on 8'^ X oo and 20^ x oo lattices. Comparison of the two volumes shows the 
features which are robust. In the broken phase the lines of constant mass 
ratio flow towards the critical point but curve away just before reaching it. 
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Figure 32: Lines of constant E^^^ jj^ (solid) and of constant Hr/ f^ (dotted) 
in the plane of the bare parameters from Monte Carlo results on an 8^ x 16 
lattice. E^^/hr = 0.7 to 1.6 in steps of 0.1 (right to left). /iR//,r = 1-5 to 
3.5 in steps of 0.5 (left to right). 



52 



m 



u. ± 


\ 1 \ 


7 ^ 


I 


y \ ' \ '/ ' / ^ ' / ' ' / 




\ \ 


/ ^ 


1 


1 / / / If / ' 




\ X 


/ 


\ / 


1 \ / / / / \ / / 





\ ^^ / 


/ 


\ / 


III V / " 






\ 


\ / 














0.08 


\ / ^ 


\ 


\\^ 1 


\ \ \ I H 1' / 














\ / \ 




\ 1 i 


j I 1 1 \ I / / ^ / 






\ 








/ \ 


\ 


/\ 


'^ I 1 \ / / / / 






\ 






0.06 


A '' 


\ 
\J 


\ I 


\\\ \// / '^ / - 




/ \ ^^ 




^J 


r 1 \ / / / \ / 




' \ 




\ A 


/ / \/ / / / / / / y/ 




\ X 


'^ 


\ / \ 


/ hi I" / / /\ 






\ 
















0.04 


y^ \ 


^^ 




\/ N A / / \^^ 




^ ^y^ \ 




^^ / 


/ \ 1 h / / y /^ /\ 




y\ y' 




\ ^/ 


/ ill! / / / y^ \ 




\ ^^ 










'^ ^^.--^ 


'^ K^ 


"'^^^ \ y^ 


\ / y / / ^J / / "" / X \ 




^^..-'^ 










^'-""'^ ^x ^.--'''^ 


^ 


^^-'-''^ 


y\ / y / /^ / / /^ y^ ^,.-'-^\ 


0.02 


^^^...--"^ 


^^..^'--' 


\---'^\ 


>( /y / X / / A/ ^^.^-^^'^ \ - 




''^^II— "-""""'"^^ 




""""^^1--^ 


^\^^ /\A / \/^X ^/A ^^^.-^"'^ ^^ _ 




— -""""^^^^^^l----''^ ^ 


J--^..^ 


— — ^S-^-^ 


"^AA^V'-^O^/^y^ y''^^^ J^<'^^^^ ^-^-"""'^^^ 


n 


^=^^ 


— ^ 


" 1 


AvAa>AAa-^^ 



0.1 



0.2 



(a) 



0.3 



0.4 



go 




0.4 



(b) 



go 



Figure 33: Lines of constant E^^ fxji (solid) and of constant jJiYi/ f-n (daslied) 
from Schwinger- Dyson calculations on 8^ x oo (a) and 20^ x oo (b) lattices. 
For E'^^ / fiYi tlie ratios are: (a) 0.2 to 1.8 and (b) 0.1 to 1.9 in steps of 0.1. 
In both figures ^k/ fir = 0.5 to 4.0 in steps of 0.5. 
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. In the 
and the 



In the symmetric phase the ratio is always just below 2 because we have 
used E^\ the lowest energy level in the vr channel. (When comparing these 
figures with fig. |32| remember that the equivalence between the Monte Carlo 
and the Schwinger-Dyson equations is not exact, there being a small shift in 
Qc between them.) 

Extrapolating to infinite volume by the Schwinger-Dyson equations gives 
fig. 0. In this limit we can identify resonances and find their masses. These 
are shown in the diagram as dotted lines. Now all the lines fiow into the 
critical point. This type of fiow diagram was also found in QED [[TT 
broken phase we can see that the lines of constant mass ratio, fig. |3^, 
lines of constant /Ur/Ztt, fig- ^ take different paths just as in fig. |3^, again 
showing that there are no lines of constant physics. 

To test for the existence of lines of constant physics in the symmetric 
phase we have looked at the ratio between itIt, and F^. These quantities are 
found from the real and imaginary parts of the pole position. In fig. |3^ lines of 
constant r.,r/mT^ are compared with curves of constant ^yi/uIt^. Once again 
the different flows cross. In this phase too there are no lines of constant 
physics. To display this information another way we show in flg. |^ the 
way in which TT^/niT^ varies along a particular curve of constant mass ratio 
{fj,R/m.„ = 0.20). Because the phase space available for the decay depends on 
the mass ratio the variation in FTr/m^ along this path is not due to kinematics, 
but must reflect a variation in the strength of the coupling between the vr 
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Figure 34: Lines of constant rriT,/ jji-^ on an oo^ lattice. Solid lines are used 
where the vr is a bound state and dotted lines where it is a resonance. The 
ratios are 0.1 to 2.0 in steps of 0.1 (bound state) and 3 to 10 in steps of 1 
(resonance) . 
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Figure 35: Lines of constant /xr/Ztt on an cxd^ lattice. The ratio runs from 
2.5 (inner curve) to 4.0 (outer curve) in steps of 0.5. 
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Figure 36: A comparison between the flows of constant T-,^/m,T^ (solid lines) 
and constant fiR/m^j^ (dotted lines) in the region where the tt is a resonance. 
f^n/fTT'TT runs from 0.1 to 0.4 in steps of 0.1, F^/ttItt from 0.1 to 0.3 in steps 
of 0.05. 
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Figure 37: The variation in T^/rriT^ along the curve /xpt/mTr = 0.20. 
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Figure 38: The variation in T^/m^ along the curve /iR/m^ = 1/2. 
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Figure 39: Lines of constant ^-^/m^ on an oo^ lattice. The ratios are 0.51 
and 0.50 (bound state and threshold: solid lines); 0.45 to 0.05 in steps of 
0.05 (resonance: dashed lines). 
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Figure 40: A comparison between the Hows of constant Tfj/mcr (solid lines) 
and constant fi-n/rrifj (dotted lines) in the region where the a is a resonance. 
f^n/iT^'a runs from 0.05 to 0.45 in steps of 0.05, Tfj/mcr from 0.0 to 0.3 in steps 
of 0.05. 
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and the fermion. As we approach the critical point, /ir small, this physical 
coupling decreases, consistent with the hypothesis that the NJL model is 
trivial. 

Note that on the curve m,r = 2yUR, the threshold for pion decay, both 
Ftt and fjr vanish. This raises the possibility that this threshold could be 
a line of constant physics. To check this we have to look at other particles 
such as the a. In fig. ^ we plot the variation of the ratio T^/m^ along the 
line m-TT = 2/ipj,. We see that the decay rate of the a is not constant along 



this curve, but decreases slowly as we approach the critical point. Fig. |3^ is 
in fact reminiscent of fig. ^ and again consistent with a trivial continuum 
limit for the NJL model. 

In fig. |3^ we show the ratio of fermion to a mass as found from the 
Schwinger-Dyson equations. Above Qc there is a large region where this ratio 
is always very close to the value of 0.5 found in 0. In the symmetric phase the 
cr is a resonance with a mass similar to that of the vr, as would be expected 
from unbroken chiral symmetry. (Because the a mass hovers around 2//r 
in the broken phase the mass ratio E^^/fi^i is essentially constant and a 



Monte Carlo flow picture shows nothing new.) In fig. EOwe compare the flow 



patterns for the ratios Tcr/m^ and jji-^/nia- As in fig. ^, the corresponding 
picture for the pion, we see a crossing of the two flows, once again a sign of 
non-renormalisability. 
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7 Conclusions 

We have made a thorough investigation of a lattice version of the NJL model 
using both the Monte Carlo method and Schwinger-Dyson equations. The 
interplay between both methods allowed us to come much further than with 
either method alone. We extensively used the Schwinger-Dyson equations 
because of rather large finite size effects. 

Since we were interested in the chiral symmetry properties we worked with 
staggered fermions. A phase transition was seen at about Qc ~ 0.280(4). In 
the Goldstone tt channel we have successfully identified a bound state in 
the broken phase and a resonance in the symmetric phase. In the a and p 
channels we have seen resonances. We have found that particles in the vr 
and a channels become massless at the critical point while particles in the p 
channel scale with the inverse lattice spacing. 

Equipped with the results from spectroscopy we worked out renormali- 
sation group flows in the bare parameter plane. We do not find any lines 
of constant physics: lines of constant Pr/ fn cross lines of constant m^//iR. 
The ratio nicr/pR has come out essentially constant in the broken phase and 
so has no bearing on this question. The absence of renormalisability implies 
that one should be cautious in applying renormalisation group techniques as 
is sometimes done in the top-mode standard model. 

The NJL model has attracted speculations such as a non-Lorentz invariant 
vacuum or a massless vector state 0. These do not appear to be realised at 
the phase transition that we investigated. 

Our lattice model can of course be embedded in a generalized NJL model 
with more bare parameters |^. It might then be possible to find regions in 
this extended parameter space where the model is (weakly) renormalisable. 
If such a region intersects the part of the parameter space we have studied 
our physical results (absence of light vector states, occurrence of resonances, 
etc.) would hold inside this whole region. 

We have now reached a comprehensive understanding of the four dimen- 
sional NJL-Model concerning the chiral condensate, the fermion mass, meson 
spectroscopy and renormalisation properties. 
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A Schwinger-Dyson Equations for the Com- 
posite States 

In this appendix we shall derive meson correlation functions using the Schwin- 
ger-Dyson equations shown in fig. O. K is the two particle irreducible kernel, 



which is in general a function of 4 coordinates, and so on a lattice of volume 
V has V^ values. This is reduced by a factor V due to translation invari- 
ance. If we consider a fermion-antifermion pair with definite centre of mass 
momentum p, momentum conservation further reduces this to V'^ values. K 
can therefore be regarded as a V^ x V^ matrix that links the separation vec- 
tor of the incoming fermion-antifermion pair (with V possible values) to the 
relative distance of the outgoing pair. This becomes much more tractable 
when we replace the full kernel by the bare four fermion interaction. The 
resulting equations are at the same level of accuracy as the order g^ Schwin- 
ger-Dyson equations which were used for the fermion propagator and the 
chiral condensate. 

Because the bare kernel is very short range, involving only fields separated 
by distance or 1, we only need the two particle wave function at 9 values 
of the separation (0 or ±/i). This allows us to write the Schwinger-Dyson 
equations in a 9 x 9 matrix form, so we can express the solution of them in 
terms of the momentum space meson propagators 
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and i]f^{x) is defined in ( p.3|) . In this basis the four fermion kernel K is given 
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and the Schwinger-Dyson equation reads 

P = M + MKP . 



(A.4) 



As can be seen from the figure M is given by the independent propagation 
of a fermion and antifermion, i. e., the t-channel disconnected part of P 



Mij{p) = — \Oi{—p)0 j [p)) t-cha.nnel disconnected • 
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The "bubbles" (see fig. 0) in M are 



(A.6) 
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Figure 41: Graphical representation of the "bubbles" in eqs. ( 1A.6|) and ( |A.7f) . 



and A^, F^ are given in (p.l4|) . Equation (|A.4|) is solved by 



P ={1-MK)-'M. (A.9) 

This is the full solution including the annihilation term. The matrix inversion 
is done numerically. 

Because we are working with staggered fermions we get mesons at rest if 
the spatial components of p are all either or tt. Other values of p give moving 
mesons and can be used to check the restoration of Lorentz invariance. 

In most of our Monte Carlo calculations we have measured only the 
fermion line connected part of the meson propagator as is customary. This 
means that t-channel exchange of auxiliary fields is kept, but s-channel (i. e., 
annihilation) auxiliary fields are dropped (see eq. ( |4.6| )). When we make this 
same approximation in the Schwinger-Dyson equations only the first term 
in the kernel is kept. We no longer have to consider bond operators, so the 
9x9 matrix becomes a 1 x 1 matrix. The meson propagator simplifies to 



P{p) ^ E(fe)Wfe)(0))e-^^-^ 

X 

(^ 

l + a(p)2fl'oE/.cosp^' 



(A. 10) 
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The propagator P{p) includes the first four multiplets shown in table |^ be- 
cause the s-factors in the table are all of the form exp{ip-x), i. e., all represent 
a simple shift in momentum in eq. ( |A.1CI| ). We are most interested in the prop- 
agators Cj- (t) for stationary mesons, i. e., mesons for which p = OmodTr. Let 
the Fourier transform of C^- (t) be C^ (cu) (see eq. (|4.13| ) for normalisation). 



The channel i = 1 (the n' /a channel) corresponds to 

C'W(cj) = P(0,0,0,7r + cj), (A.ll) 

i = 2 (the TT channel) corresponds to 

C'f^(cj) = P(7r,7r,7r,7r + cj), (A.12) 

and the p'/a and p/b channels {i = 3, 4) to 

C^iuj) = P(7r, 0,0,71 + uj) = P{0, 71,0,71 + to) = P{0, 0, TT, vr + u) (A.13) 
and 

Cf )(CJ) = P(0, 71,71,71 + UJ) = P{7X, 0,n,7l + Uj)= P{7X, TT, 0, TT + Uj) (A.14) 

respectively. Representing the C"s as particular values of P{p) shows that 
the fact that the different i and k channels don't mix is simply a consequence 
of the conservation of 3-momentum. Inserting the above momentum values 
in ( |A.10| ) gives the results stated in eq. ( [4.14|) . 
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B The Infinite Volume Limit of the Meson 
Propagators 

In this appendix we study the infinite volume hmit of the meson propagators 
( ^■141 ). We have aheady taken the Lj — »• cxo hmit in ( [4.20[ ). We now let 
Ls — > oo. In the infinite volume limit the different boundary conditions in 
the space and time directions no longer matter and so Fg = Ft = F. 
The fermion loop integrals in ( [4.18|) are all of the form 



"^"^ dki dk2 dky, ^^ f . 2, .2; -2 



27r 27r 27r 



■/i(sin ki + sin k2 + sin k^ 



(B.l) 



In the continuum we would naturally simplify such a 'spherically symmetric' 
integral by changing variables to k\ + k^ + k'^. In the same way we use the 
variable 

r = -3 + 2(sin2 k^ + sin^ ^2 + sin^ k^) (B.2) 

to simplify the lattice integrals. (The shift and normalisation make some 
of the following expressions simpler than they would otherwise be.) The 
integral (p.l|) then becomes 



-3 



drhi — I — r 



(B.3) 



The 'phase space' function is 

r2^ d^k 



(27r)3 
+°^ du 



5(r + 3 - 2(sin^ fci + sin^ ^2 + sin^ k^)) 



Ztx 



(B.4) 



is outside the range (—3,3) and has branch points at r = ±1. These 
branch points are relevant to the analytic structure of the spectral density. 

Expressing the A's and 5's in eqs. ( ^.15| ) and ( |4.1(j| ) in terms of r and 
using eq. (|2.15| ) to eliminate N gives 



A^ -- 


= F2(sinhVR + 


{r + 3)/2) , 


Bl - 


= F^l^sinhVR- 


(r + 3)/2) , 


Bl - 


= F2(sinhVR + 


(r + 3)/2) , 


Bl - 


= F2(sinhVR- 


(r + 3)/6) , 


Bl - 


= F'(sinhVR + 


(r + 3)/6) . 



(B.5) 
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We have now reduced our original four- dimensional integrals to one dimen- 
sional integrals from which we can find the analytic structure of ai{uj) and 
therefore also of Ci{u!). In an infinite volume ai{uj) is given by an integral of 
the form 

a.M = Trj'dr<p{r)(-—^J±^ 

r^J-3 \2smh /ir + r + 4 — coscu 

^J^ ] . (B.6) 

2 sinh /iR + r + 4 + cosujj 

f+{r) and f-{r) are analytic (and real) in (—3,3). The ensuing singularity 
structure is sketched in fig. |l^ and discussed in the main text. 
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